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The Kitaev-Heisenberg (KH) model has been proposed to capture magnetic interactions in iridate 
Mott insulators on the honeycomb lattice. We show that analogous interactions arise in many 
other geometries built from edge-sharing IrOe octahedra, including the pyrochlore and hyperkagome 
lattices relevant to Ir2 04 and Na4lr308 respectively. The Kitaev spin liquid exact solution does not 
generalize to these lattices. However, a different exactly soluble point of the honeycomb lattice 
KH model, obtained by a four-sublattice transformation to a ferromagnet, generalizes to all these 
lattices. A Klein four-group = Z2 x Z2 structure is associated with this mapping (hence Klein 
duality). A finite lattice admits the duality if a simple geometrical condition is met. This duality 
predicts fluctuation free ordered states on these different 2D and 3D lattices, which are analogs 
of the honeycomb lattice KH stripy order. This result is used in conjunction with a semiclassical 
Luttinger-Tisza approximation to obtain phase diagrams for KH models on the different lattices. 
We also discuss a Majorana fermion based mean field theory at the Kitaev point, which is exact on 
the honeycomb lattice, for the KH models on the different lattices. We attribute the rich behavior of 
these models to the interplay of geometric frustration and frustration induced by spin-orbit coupling. 



I. INTRODUCTION 

Long viewed as a perturbative correction, relativistic 
spin-orbit coupling has in recent years been increasingly 
asserting its role within condensed matter physics. It 
took center stage with topological insulators, time rever- 
sal invariant states of electrons with no strong interac- 
tions that use spin-orbit coupling to generate nontrivial 
topology in the band structure^l'Sl. Electron correlations 
may amplifjBUl the effects of spin-orbit coupling (SOC), 
enriching the taxonomy of possible phases. Thus Mott in- 
sulating states of heavy magnetic ions could realize novel 
Hamiltonians, which may be hitherto unexplored or not 
thought to describe real materials. 

One such S=l/2 Hamiltonian has been proposed by 
Jackelli and KhaliuUin^ to occur in the honeycomb iri- 
dates Na2lr03 and Li2lr03. It includes the Kitaev ex- 
change, a nearest neighbor ising coupling of spin compo- 
nent 7 € {x, y, z} set by the spatial orientation of the 
bond^ ^. The pure Kitaev honeycomb Hamiltonian is ex- 
actly solvable with a quantum spin liquid (QSL) ground 
state of a gapless Majorana coupled to Z2 fluxegS). The 
proposed magnetic model for these iridates includes the 
Kitaev as well as SU(2) symmetric Heisenberg coupling, 
yielding the Kitaev-Heisenberg S — 1/2 Hamiltonian . 
We parametrize it as 

hkb = Y.'i [(1 - • - 2«<57-^;-] (1) 

witlPl?? = ±1 and -1 < a < 1. Here 77 sets the sign of 
Heisenberg exchange, and negative a involves the same 
sign for both exchanges. Pairs of endpoints of the two 
a-segmcnts for rj = +1,-1 are identified as a single point 
by the product rja — +1 (FM Kitaev) and similarly rja = 
— I (AF Kitaev), forming an (?7, a) parameter ring. We 



will primarily focus on this idealized Hamiltonian but 
also consider some extensions such as farther neighbor 
couplings. 

The phase diagram of ifxH on the honeycomb lat tice 
is known from a combination of exact diagonalizatiorF'Sl 
and the presence of exactly soluble points. In addition 
to the exact solution using majorana fermions of the 
Kitaev Hamiltonians a = ±1, and the obvious SU(2)- 
symmetric ferromagnet {j] — — l,a — 0), a four sub- 
lattice site-dependent spin rotatiort^ transforms Hkh at 
77 = -1-1, a = 1/2 into a ferromagnet in the rotated basis. 
The original spins are then "stripy" ordered. Neel order 
from the Heisenberg antiferromagnet is unfrustrated on 
the bipartite honeycomb, and was recently shownii* to 
map under this transformation to a physical parameter 
regime hosting the spin pattern known as "zigzag" . 

This zi gzag phase was determined in recent 
experimentalSHni to be the low temperature order- 
ing pattern of Na2lr03. The zigzag order was earlier 
theoretically found to be most stabilized by combining 
Kitaev interactions and the further neighbor exchanges 
J2 , J3 whi ch naturally arise across a honeycomb 
hexagorP^E^, and which together fit the available exper- 
imental data in comparisons to exact diagonalization'^l. 
Indeed, within a classical approximation to the phase 
diagram (Luttinger-Tisza described below), the zigzag 
phase within the pure Kitaev-Heisenberg model lies 
nearly at the boundary of the large zigzag-ordered 
region stabilized by J2 , J3 exchanges^. Interestingly, 
the zigzag phase in its Ji — J2 — J3 limit and in its Eq. [T] 
limit may o ffer ex perimentally relevant distinguishing 
characteristici^^^. 

So far only the honeycomb iridates Na2lr03 and 
Li2lr03 have been studied in the context of the Kitaev- 
Heisenberg model. Despite initial worries that trigonal 
distortion would invalidate the derivation of the Kitaev 
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exchange discussed below, recent resonant inelastic x-ray 
scattering resultJ^ support the validity of the strong 
spin-orbit coupling approach. For the sodium iridate 
Na2lr03, attempts to extract the magnetic Hamiltonian 
from fits to experiments i ncluding susceptibility and spin 
wave spectrgfliffi^IIIIlIIlll qj^^ electronic propertieJ^^HSl 

have so far proved unable to distinguish between sub- 
stantial Kitaev exchange and a complete lack of it. Few 
experimental results on magnetic behavior in the lithium 
iridate Li2lr03 are currently available, though the rela- 
tively small magnitude of the Curie Weiss scale extracted 
from susceptibility suggests the Kitaev exchange may be 
stron^i^^ . 

Beyond the possible Kitaev-Heisenberg physics in the 
layered honeycomb iridates, other iridates have also at- 
tracted much attention. Layer ed compounds include 
the M ott in sulator Sr2lrOj2il25| and its bilayer variant 
SrglrOTiaMZI, both with ordered SOC magnetic moments. 
Notable examples with a fully three dimensional struc- 
ture include the 2-2-7 pyrochlore iridates, where chang- 
ing the A sit e rar e earth metal yields radically varying 
properties* ^ the sodium iridate Na4lr308 spin liq- 
uid candidate, with Ir on the pyrochlore-descendent hy- 
perkagome lattic^^oJ and a recently epitaxi ally-stabilized 
iridium spinel Ir2 04 with empty cation sites^D leaving Ir 
on a pyrochlore lattice. Despite the variety of elemental 
composition and geometrical structure in this list, there 
is a simple but fundamental distinction separating the 
latter two compounds from the others listed. 

In this manuscript we show that the iridates Na4lr3 0g 
and Ir2 04, as well as possible compounds in certain other 
geometries, may be described by generalizations of the 
Hamiltonian Eq. [l] to the relevant lattices (hyperkagome 
for Na4lr308 and pyrochlore for Ir204). The key quan- 
tum chemistry ingredients which can generate the inter- 
actions -/Jkh have been already pointed out by Jackeli 
and KhaliuUirP but the extension to three dimensional 
lattices, as well as to these compounds, has not been 
previously exposed. We begin by recalling the derivation 
of i?KH and systematically extending it to other geome- 
tries in two and three dimensions; it applies when oxygen 
octahedra are edge-sharing, yielding lattices that are in 
a certain sense subsets of the fee. We then proceed to 
investigate the phase diagram of -Skh on these lattices, 
using primarily analytical approaches. We generalize the 
honeycomb four-sublattice transformation into a dual- 
ity on graphs and lattices with Kitaev labeled bonds in 
any dimension; we shall refer to it as the Klein dual- 
ity since, as we shall show, it is structured by the Klein 
four-group = Z2 x Z2. We give a simple algorithm de- 
termining which graphs admit the duality, based on this 
Klein group structure. The Klein duality gives stripy 
phases as FM-duals. Diagonalizing the classical version 
of with spins of unconstrained length (i.e. the 

Luttinger-Tisza approximation), we identify unconven- 
tional ordering patterns and also find hints of quantum 
magnetically disordered phases, most interestingly on the 
hyperkagome. The Luttinger-Tisza phase diagrams are 



shown in Fig. [T] To directly capture Majorana fermion 
quantum spin liquids analogous to the Kitaev honeycomb 
QSL, we decompose spins into Majorana combinations of 
Schwinger fermions, a mean field treatment which is ex- 
act for the Kitaev honeycomb model, finding on all other 
lattices fermionic QSLs which break time reversal and 
carry gapless excitations, but which are not exact solu- 
tions. For the honeycomb and hyperkagome pure Kitaev 
Hamiltonians, the lattice fragments under a bond type 7 
into disjointed localized clusters, giving flat bands in the 
Majorana mean field as well as in the Luttinger-Tisza ap- 
proximation, which hints at a possible analogy between 
the honeycomb Kitaev QSL and the Kitaev Hamiltonian 
ground state on the hyperkagome. 

We focus on two candidate materials, while also consid- 
ering other related compounds. The recently epitaxially 
fabricated Ir2 04 is a spinel without the A cation, leav- 
ing the iridium ions on a pyrochlore lattice with oxygens 
positioned appropriately, as described below; Ir204 was 
found to be a narrow gap insulatoil^. The spin liquid 
candidate Na4lr3 08 is an iridate with S = 1/2 moments 
on the three dimensional hyperkagome lattice, which ex- 
hibits no magnetic order down to at least 2KP. In ad- 
dition to these two iridates, this study may also capture 
compounds in which iridium is replaced by a transition 
metal ion with strong spin orbit coupling, intermediate 
correlations and valency appropriate for a magnetic ef- 
fective spin one half model (see below). Recently, the os- 
mium oxides CaOs204 and SrOs2 04 were computation- 
ally predictecP^ to be stabilized in the spinel structure 
relevant to Kitaev-Heisenberg physics; if they indeed ex- 
ist in this geometry, i?KH should form at least part of 
their magnetic Hamiltonian. Kagome and triangular lat- 
tice iridates may be seen as appropriate layers within epi- 
taxially stabilized Ir2 04, coupled together in a nontrivial 
manner. Triangular lattice iridates could potentially also 
be stabilized as analogues of the cobaltate^^ Naa;Co02, 
where Co is on a triangular lattice; the preferred valency 
would exist uniformly only in the limit x — ^ 0, but small 
X should offer interesting perturbations as well as likely 
separate layers of triangular lattices. However, no trian- 
gular lattice iridate with the relevant edge sharing octa- 
hedra structure is currently available; a compound of the 
type Naa;Ir02 may or may not turn out to be stabilized. 

In addition to the honeycomb Kitaev-Heisenberg 
model, other previous work has investigated Hamiltoni- 
ans related to Hkh on other lattices. Chen and Balents^fil 
studied spin Hamiltonians on the hyperkagome lattice for 
Na4lr3 08, in the strong and weak SOC limits (relative to 
octahedral distortions). Within the strong SOC case they 
considered the single superexchange pathway via oxygen 
ions generating the single point ffKH at = 1/2, for 
which they found that classical configurations of stripy 
patterns were completely unfrustratecP^. Superexchange 
via oxygen ions generating anisotropic spin interactions 
for the SOC Kramers doublet in Na4lr308 was also con- 
sidered by Micklitz and NormarPsM! in electronic struc- 
ture computations and associated microscopic tight bind- 
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FIG. 1: Phase diagrams using the Klein duahty as well as the Luttinger-Tisza approximation (LTA). 

Phase diagrams are shown for the Kitaev-Heisenberg Hamiltonian, Hkh = V ~ \<^\) Si ■ Sj — 2riaS^^' Sj^' with 
r] = ±1, on the various lattices. Note that the parameter space is a ring: for the a parameter segments shown here 
the endpoints should be identified, i.e. writing the parameter as a), identify the points (-1,-1)?^(+1,+1) and also 
identify the points (-1,+1)=(+1,-1). Rich phase diagrams are found. 2D and 3D stripy magnetic phases (blue), 
found on all lattices, are exact and fluctuation free by the Klein duality at the FM SU(2) point (yellow star) 
77=+l,a=l/2. On the kagome and hyperkagome lattices, outside the FM SU(2) points, the FM and stripy orders 
are given non- unit-length spins by the LTA (gray shading), suggesting frustration. Extensive degeneracy of LTA 
ordering wavevectors hints at a non-spin-ordered "quantum" phase, labelled "Q ": wh ere the Hamiltonians hosting Q 
phases have been solved, exactlj^for the honeycomb at q;=±1 and numer icalljEUm for the ka gome at 77=+l,a=0, 

they have turned out to host non-magnetic phases. The hyperkagome hosts Q points at the Heisenberg 
antiferromagnet and its Klein dual, as well as at the pure Kitaev Hamiltonians, because any single Kitaev bond type 

fragments the hyperkagome into disjoint clusters (Fig. |4|. The 120° triangular lattice and Neel and zigzag 
honeycomb orders are found with normalized spins. Apparent ordering with LTA non-normalized spins is found at 
incommensurate wavevectors on the triangular lattice and in the clustcr-ferrimagnet and cluster-antiferromagnet 
(AF) regimes on the hyperkagome. Kagome, fee and pyrochlore lattices also host frustrated regimes with no 
definitive ordering within the LTA, as described in the text (gray). 



ing parametrization. Very recently the classical i?KH 
Hamiltonian was studied on the triangular lattice by 
Rousochatzakis, Rossler, Brink and Daghofei^^. Their 
study included a classical Monte Carlo computation sug- 
gesting that Kitaev exchange may stabilize an incommen- 
surate vortex lattice of the Z2 topological defects known 
to exist in the Heisenberg antiferromagnet 120° order. 



II. KITAEV COUPLINGS IN LATTICES 
BEYOND THE HONEYCOMB 

Generating the Kitaev coupling requires a subtle recipe 
with ingredients from chemistry, geom etry, a nd a hierar- 
chy of energy scales, as we now recalW^^. Spin orbit 
coupling is key, together with (intermediate) correlations; 
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let us focus on iridium. The iridium ions should retain 
their 5d electrons in localized orbitals, and exist in the 
4+ valence. Each iridium should be surrounded by six 
oxygen ions (or other electronegative ions with valence p- 
orbitals), which form the vertices of an octahedron cage, 
shown in Fig. [2j The octahedral crystal field splits the 
5d orbitals into an empty Cg pair and a triplet of or- 
bitals with five electrons and one hole. Strong spin-orbit 
coupling further splits e2g down to a half-filled Kramer's 
doublet, the spin-1/2 degree of freedom defining the low 
energy manifold. The final key ingredient is the geomet- 
rical structure: edge sharing octahedra with 90° Ir-O-Ir 
bond angles. 

In perturbation theory from the Mott insulator limit, 
virtual hopping of holes from iridium orbitals through 
intermediate oxygen p orbitals generate the low energy 
spin Ham iltonian. There are multiple relevant exchange 
pathJ^EIl^ When holes hop through intermediate oxy- 
gens and meet on an iridium d-orbital, the resulting 
coupling is a pure Kitaev term, and is proportional 
to JhI {{Ud - iJnWd - Jh)) « Jh/UI The iridium 
Coulomb exchange Ud and Hund's rule coupling to- 
gether specify all of the multi-band interaction parame- 
ters, due to the symmetries of d-orbitals. A second ex- 
change path, with two holes meeting on an oxygen or 
cycling around the Ir-0 square, contributes a combina- 
tion of Kitaev and Heisenberg couplings equal to i?KH at 
a — 1/2, with a coefficient and sign depending on the 
oxygen p-orbitals charge-transfer gap and Coulomb re- 
pulsion. Direct iridium wavefunction overlap gives a pure 
Heisenberg coupling. RecentljEl^ an additional pathway 
through the higher eg orbitals has been proposed to be 
relevant as well, contributing i?KH at 77 = —1, a — 1/2. 
The interplay of these exchanges suggests a may not be 
computable microscopically. 

Generalizing this derivation to geometries beyond the 
honeycomb requires preserving the edge sharing octahe- 
dra with 90° Ir-O-Ir bonds. Many commonly studied 
iridates such as the layered perovskites and the "2-2- 
7" pyrochlores have corner sharing octahedra and thus 
are not captured by this derivation. Fig. [2] shows two 
adjacent octahedra, with edges color coded by the spin 
component coupling they generate when the octahedra of 
neighboring Ir ions share that edge. It is evident that all 
twelve octahedra edges may be shared while still main- 
taining 90° bonds and three-fold symmetries (coupled 
space and spin rotations). Tiling octahedra which touch 
along edges builds a face-centered cubic (fee) lattice of 
the octahedra centers. 

We thus find that in two and three dimensions, all lat- 
tices whose graph of nearest neighbor bonds is a subset 
of the nearest neighbor bonds of the fee, including the 
fee itself, may host analogues of the Kitaev exchange. 
Possible geometries include the kagome and triangular 
lattices in two dimensions, and the face-centered cubic, 
pyrochlore (as realized in spinel-based compounds) and 
hyperkagome geometries in three dimensions. These are 
shown in Fig. [Sj The materials discussed above all have 




FIG. 2: Edge-sharing IrOg Octahedra generating 
Kitaev exchange. Iridium ions (spheres) are each 
coordinated by six oxygen ions forming vertices of 
octahedra. Octahedra of neighboring Ir ions share 
edges. Dotted (purple) lines show the 
iridium-oxygen-iridium hopping paths, which form a 
square with 90° angles. As described in the text, these 

superexchange paths generate an ising interaction 
between the iridium effective spins, which couple a spin 
component x, y ov z depending on the orientation of the 
shared octahedra edge (shown in red, green and blue). 
Ir lattices hosting this Kitaev exchange must arise from 
a regular tiling of these edge-sharing octahedra. 



these edge-sharing octahedral structures and their mag- 
netic Ir ions form one of these lattices. As for the honey- 
comb iridates, reduced crystal symmetry distorting Ir-0 
angles away from 90° will generate other magnetic ex- 
changes. Despite apparently strong > 10% distortions of 
the bond angle in the sodium honeycomb iridate and a 
slew of experiments on this material, a Kitaev exchange 
comparable to or even stronger than the Heisenberg ex- 
change is still consistent with current experimental re- 
sults, suggesting a hopeful outlook for the other materi- 
als. 

Note that the quantum chemistry considerations pic- 
tured in Fig. [2] tightly constrain the possible lattice re- 
alizations of Eq. [1] Specifically, these constraints are 
tighter than those imposed by naive symmetry consid- 
erations of SOC. For example, it is natural to define an 
implementation of SOC that couples spin component 
to bonds along z, i.e. locks the Bloch sphere to real 
space. This would generate Eq. [T] on the simple cubic 
lattice with S^S^' coupling along z bonds, as well as on 
the square lattice with 7 = a;, y. But the exchange path- 
ways of Ir t2g orbitals forbid this scenario. Instead, the 
analysis above shows that for t2g orbitals as in iridium, 
SOC couples spin component to the fee lattice bonds 
lying normal to z. The simple cubic lattice version of 
iJKH cannot be generated, and a compound structured 
as layers of a square lattice would collapse its Kitaev ex- 
change to uniform ising couplings along all square lattice 
bonds. 

The honeycomb and hyperkagome lattices share a com- 
mon feature distinguishing them from the other lattices: 
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FIG. 3: Kitaev-Heisenberg lattices. Iridium ions 
arranged in these lattices may generate the Kitaev spin 
exchange, coupling component x, y, z on bonds colored 

red, green, blue respectively. A blue colored bond 
connecting two Ir sites implies that the respective IrO^ 
octahedra share a blue edge as in Fig. [2j We also list 
examples of possible relevant iridium compounds which 
form these lattices. 



if we only keep bonds of a single Kitaev type 7, the lattice 
fragments into localized disconnected clusters. On the 
honeycomb, each cluster contains two sites, and forms 
the unit cell. On the hyperkagome, each cluster contains 
three sites, arranged into a line segment. For a given 
bond label 7, the twelve site unit cell fragments into four 
disjointed clusters, whose line segments are oriented par- 
allel within each of two pairs and perpendicular between 
the pairs. The structure on the hyperkagome unit cell is 
shown in Fig. |4] As discussed below, this fact has dra- 
matic repercussions for the Kitaev Hamiltonians in both 
the Luttinger-Tisza approximation and in the Schwinger 
fermion Major ana mean field (which on the honeycomb 
describes the Kitaev QSL). In both cases, certain excita- 
tions only propagate within a single Kitaev bond type 7, 
and the localized disconnected clusters imply these exci- 
tations must have completely flat bands. 




FIG. 4: Hyperkagome unit cell and 
decomposition into Kitaev clusters. A symmetric 

depiction of the hyperkagome structure, highlighting 
the four disjointed three-site clusters which split up the 
unit cell for a given Kitaev bond label (red,green,blue). 
The fact that the lattice fragments into these disjointed 

clusters for a given Kitaev bond type has substantial 
repercussions, as described in the text. The four clusters 
of a given type appear in two parallel-orientation pairs 
which are perpendicular to each other (with different 
shading). The unit cell is composed of the 12 sites 
participating in the four drawn triangular faces, as well 
as all 24 drawn bonds. We label the 12 sites by a letter 
A,B,C as in Ref . fiT] so that A spins lie on midpoints of 
type-A (here blue) clusters, etc; and by a number 1-4 
(chosen to not repeat within a triangle face or cluster). 

The camera angle, i.e. the vector pointing into the 
page, is just slightly off (up-right) of a Cartesian (here 
also Bravais) axis; for Na4lr308 it is the vector from an 
iridium ion to a neighboring coordinating oxygen. 



III. KLEIN DUALITY AND HIDDEN 
FERROMAGNETS 

A. Connections to previous work 

Exactly solvable quantum Hamiltonians are rare in di- 
mension higher than one. It is quite remarkable that the 
stripy phase at = a = 1/2 found for the honey- 
comb Kitaev-Heisenberg modeP is exact, a hidden fer- 
romagnet exposed by the site dependent spin rotation 
which quadruples the unit celP^. Unlike Neel order on 
even bipartite lattices, this stripy antiferromagnetic or- 
der is exact and fluctuation- free at a = 1/2. 

This "four sublattice rotation trick" has been known 



by Khaliullin and Okamoto for t2g orbitals in a cu- 
bic environment since as early as 200^^. It was used 
for Kitaev- Heisenberg-like Hamiltonians in ferromagnetic 
titanateS^^iMi as well as in other systems, incl uding an ex- 
plicit transformation on the triangular lattice^^ESl to find 
the dual of 120° order for C0O2. It was then applied to 
the honeycomb lattice by Chaloupka, Jackeli and Khali- 
ullin in their derivation of the Kitaev-Heisenberg model 
for the honeycomb iridates^. However, its general struc- 
ture has not been previously elucidated. We will now 
show that this duality transformation may be defined 
on general graphs with Kitaev 7 bond labels and that 
it has the structure of the Klein four-group, isomorphic 
to Z2 X Z2. This will then lead to a geometrical con- 
dition specifying which lattices and finite graphs admit 
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the Klein duality, a result especially useful for designing 
finite graphs for numerical studies. 

B. Deriving the Klein duality on Kitaev lattices 

We begin by defining a general unitary transforma- 
tion, and then we will show that under certain special 
conditions it acts as a duality transformation on Eq. [T] 
Throughout this paper, by a "duality transformation" we 
refer to a mapping between Hamiltonians that maps a set 
of Hamiltonians (and the associated phase diagram) to 
itself (of course without mapping each particular Hamil- 
tonian to itself). Consider a lattice or finite graph in any 
dimension which connects S = 1/2 spins, and assume 
each bond carries a Kitaev type label 

e {1, x, y, z}. (2) 

The set 7 G {x, y, z} corresponds to Kitaev coupling 
S^S'J on that bond, where {x,y,z} identifies a set of 
orthogonal axes in the spin Bloch sphere. The Hamilto- 
nian on the bond may have other terms such as Heisen- 
berg coupling and various anisotropies; but the transfor- 
mation will turn out to be most useful if the coupling 
includes only Kitaev and Heisenberg terms, as in Eq. [T] 
The label jij — 1 can be assigned to a bond that does 
not have a Kitaev exchange (only Heisenberg and pos- 
sible anisotropies), such as a second or third neighbor 
interaction. In general such farther neighbor interactions 
supplementing 7?kh will frustrate the transformation, so 
when making use of the Klein duality the lattice should 
usually be considered to be just the pure nearest neighbor 
Kitaev-Heisenberg model TJkh , where all bonds carry 
7 G {x, y, z}. But we will show below that certain farther 
neighbor Heisenberg interactions do preserve the duality 
structure, and may be fruitfully included as 7 = 1. 

Let us proceed by describing the relevant transforma- 
tions on individual sites. Assign each site a label 

a, £ {1, X, r, Z} (3) 

which will specify a unitary transformation on that site, 
specifically rotation by tt around the Bloch sphere axis 
5° for a E {X,Y,Z}, and no rotation for the identity 
element a = 1. Note that tt rotation around 5° flips the 
sign of the spin components perpendicular to a, so that 
the rotation ai = Z multiplies the (x, y, z) components 
of Si by the sign structure g[Z] ~ (—1,-1,1), and also 
that g[l] = (1, 1, 1). 

Now, observe that both bond labels jij and site la- 
bels Qi may be interpreted as elements of the single set 
{1, X, Y, Z}. We may turn this set into a group by defin- 
ing a multiplication rule. A possible definition is sug- 
gested by the multiplication of the associated sign struc- 
tures g, which entails for example g[Ar](;[y] = g[Z], sug- 
gesting we should define X Y = Z. The resulting multi- 
plication table is defined by 

X^ =Y^ = Z^ = XYZ = 1 (4) 



with 1 acting as the identity. This is the presenta- 
tion of the group with generators (X, Y, Z) and rela- 
tions {X'^, Y^, Z^,XYZ), known as the Klein four- group. 
The Klein group is abelian and with four elements is the 
smallest non-cyclic group; it is isomorphic to Z2 x Z2. 

There is an alternative, geometrical, way to define mul- 
tiplication on the elements and 7^^. We define the 
geometric multiplication (*) of a site i and one of its 
bonds (i, j) to be the site reached by traversing the bond, 
i* (i, j) = j. The associated Klein group elements Oi and 
"fij inherit this geometric multiplication as 

a-i * 7i,j = flj- (5) 

The Klein group product ( x ) and the geometric multipli- 
cation (*) are consistent on a bond if they give the same 
answer, ai * "fij = ai x 7^^ . We say the transformation 
given by site labels {at} is the Klein transformation if the 
geometrical multiplication is consistent with Klein group 
multiplication on every bond in the lattice. 

If the transformation site labels a^, Oj across a bond 
are consistent with the Klein group product, i.e. 

a-i X 7ij = aj (6) 

or equivalently (since elements in the Klein group square 
to the identity) 

Oi X aj = 7,;j , (7) 

then the transformation changes the form of a Kitaev- 
Heisenberg coupling in an especially simple way. This is 
simply because the sign flips g[a] multiply by the Klein 
group rules, so the diagonal spin exchange J2a ^Lj'^i^f 
transforms by 

J"j 5a[a,]<7aK] J"j- = gh.j]J"y (8) 

The components of J perpendicular to the bond type 
label flip sign, i.e. the Heisenberg coefficient flips sign and 
the Kitaev coefficient gains twice the (old) Heisenberg 
coefficient. 

Even if the transformation labels on two sites are con- 
sistent with the Klein group product on that bond, it 
might seem improbable that the Oi rotation labels can be 
chosen across the entire lattice in a pattern that is Klein 
group consistent on all bonds. Such consistency for all 
bonds is necessary for the transformation to change the 
Hamiltonian uniformly. Now the Klein group structure 
shows its worth. The condition on the transformation 
{oi} — consistency between geometric and Klein group 
multiplication on each bond — can be expressed as a 
condition which refers only to the lattice: that the 7^^ 
encountered in any closed path multiply to the identity 
1. In other words, all closed loops on the lattice must be 
composed of the identity operators 1, X"^ ,Y'^ , Z"^ , XY Z . 
Then the transformation may be consistently defined by 
Klein group multiplication of bond labels on a any path. 
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C. Geometrical condition for the Klein duality 

We have shown that the existence of the Klein duality 
can be expressed as a condition on the lattice. Using 
the Klein group structure, we can write this condition 
as follows. Any closed loop, containing x-honds, Ny 
y-bonds and z-honds, must satisfy 



Nx, Ny, Nz all even or all odd. 



(10) 



The three NiS can be all even because Klein group ele- 
ments square to the identity, or all odd because XYZ = 
1. If this condition Eq. [lO]is satisfied on all closed loops 
then the Klein duality can be constructed consistently as 
follows: choosing a reference site i which for simplicity 
will be unchanged in the duality, Oi = 1, assign any site 
j a rotation label aj as simply the Klein group product 
of the Kitaev bond labels 7 on any path from i to j. 



The constraint Eq. 10 ensures this duality construction 



is consistent regardless of the choice of paths i to j. The 
Klein duality then maps Eq. [T] to itself, transforming the 
parameters a and 77 according to Fig. [5] 

It is easy to see that the Klein duality indeed exists 
on all of the infinite lattices shown in Fig. [Sj because the 
Klein group is abelian, it is sufficient to check that the 
condition is satisfied on small local loops. For example, 
triangle faces have = Ny = N^ = 1. The condi- 
tion also holds on other lattices such as the simple cu- 
bic that can host symmetric Kitaev exchange but cannot 
generate it via t2g-p orbital superexchange. Adding pure 
Heisenberg (7 = 1) further neighbor interactions gener- 
ally spoils the Klein duality, though if all the resulting 
loops satisfy Eq. [TOj the Klein duality survives unscathed 
and moreover does not modify the pure Heisenberg 7 = 1 
interactions, even while it flips the sign of Heisenberg in- 
teractions on Kitaev-labeled bonds. This occurs, for ex- 
ample, with J3 Heisenberg exchanges on the honeycomb 
lattice, connecting sites on opposite corners of a hexagon. 

For graphs and finite lattices with periodic boundary 
conditions (PBC), considering small local loops is insuf- 
ficient; loops traversing the PBC may break Eq. [TO] and 
spoil the Kitaev duality. Such winding loops must be 
checked explicitly. Here the condition Eq.[lO]should serve 
much practical use, as finite sized versions of the Fig. [3] 
lattices with PBC are useful for numerical studies, and 
it may otherwise be difficult to construct or identify the 
choice of PBC which admit the Klein duality. In many 
cases the appropriate PBC involve nontrivial twists that, 
in a continuum limit, appear as cutting and gluing oper- 
ations on the boundaries. 



D. The Klein duality on J/kh, Klein self-dual 
points and Klein Z2 symmetry 

In order to describe the action of the Klein duality on 
the parameter space of Eq.[T], let us first discuss our 77, a 
parametrization. In Eq. [T] , the sign ry = ±1 is the sign 



of the Heisenberg exchange and the sign of a is minus 
the relative sign between the Kitaev and Heisenberg ex- 
changes. Restricting to n — +1, < a < 1 gives the 
original parameter spac^' with antiferromagnetic (AF) 
Heisenberg and ferromagnetic (FM) Kitaev interactions. 
It is clear that both the FM and AF pure Kitaev Hamilto- 
nians are each described by two parameter points, which 
must be identified. 



(77 = 



fl,a 
-1, a = 



(77 

iv = 



-1, a 

-l,a-- 



-1) 



(11) 



Identifying (i.e. gluing) these pairs makes the (77, a) pa- 
rameter space into a circle. In the axis shown at the top 
of Fig. [sj the two a segments (for 77 = -|-1, —1) are con- 
nected both in the middle where they are drawn to almost 
touch and also at their distant endpoints (where arrows 
are drawn). Comparing to the angular parametrization 
presented in Ref. [TTJ 77 = +1(77 = —1) is the right (left) 
side of the circle, and a = —1, -|-1, — 1, -|-1 increases 
going clockwise. We will also sometimes refer to both 
FM and AF pure Kitaev Hamiltonians simultaneously, 
in which case the notation a = ±1 is unambiguous. 

Now we may discuss how the Klein duality acts on 
Eq. [1] . In other words, the Hamiltonian Eq. [T] with 
certain parameters 77, a is equivalent to the Hamiltonian 
Eq. [1] on the rotated spins but with different parameters 
rj',a'. The duality is shown by the blue lines in Fig. [5] 
Note that where the blue lines are roughly vertical, the 
duality approximately just fiips the sign of both rj and 
a, i.e. just fiips the sign of the Heisenberg term. In 
general it flips the sign of the Heisenberg term but also 
adds twice the (old) Heisenberg term to the Kitaev term: 
J'jj = -Jh, J'k = Jk + 2Jh. Note that with Eq. [l]as 
written, changing a also changes the overall energy scale; 
this can be avoided by dividing Eq.[l]by (1 — |q;|), so that 
the magnitude of the Heisenberg term remains fixed at 
1. 

The duality always takes 77 — ^ —77, but acts on a in 
a nonlinear way, approximately shown by the changing 
slope of the blue lines in Fig. [5] The relation between a 
and a' is given implicitly by 



(1 - Cq')(1 - Ca) , 



a<0 

2a a>0 • ^ ' 



As is clear from Fig. [5] the relation may be written ex- 
plicitly as a piecewise function, 



-1 < a < 

l~2a n ^ „ ^ 1 

i - 2 i < a < 1 

a 2 — — 



(13) 



The prototypical duality is a strong-coupling to weak- 
coupling duality, such as the duality relating the para- 
magnetic and ferromagnetic phases in the transverse 
field (quantum) ising model. The Klein duality is not 
amenable to this interpretation. First, there is no sense 
of weak and strong coupling regimes within the parame- 
ter space of Eq. [T] . Second, this parameter space forms 



8 



I A7 = +1: AF-H 1 I A7 = -1: FM-H 1 

a -1 -0.5 0.5 1 -1 -0.5 0.5 1 




1 0.5 -0.5 -1 1 0.5 -0.5 -1 a 

I A7 = -1: FM-H 1 I A7 = +1: AF-H 1 

■A', SU(2) symmetry 



FIG. 5: Action of the Klein duality on the Kitaev-Heisenberg Hamiltonian. The ij, a parameters on the 
top line (oriented left to right) of the Hamiltonian Eq. [Tjmap according to the blue lines to 77, a parameters on the 
bottom line (oriented right to left) in the Klein-dual Hamiltonian. Thick blue lines map the points shown exactly, 
thin blue lines are a qualitative sketch. Note that the right and left edges of the figure are identified, forming a ring. 
Both pure Kitaev Hamiltonians (a = ±1) are self-dual, mapping to themselves. The points at a = 1/2 and rj = 
rj = —1 are dual to the SU(2) symmetric FM-Heisenberg and AF-Heisenberg points (yellow and brown stars) 
respectively, with the FM dual point hosting the exactly soluble stripy phases on all lattices in Fig. [31 



a ring, and rather than a single self-dual point, it offers 
two distinct Hamiltonians which are self-dual under the 
Klein duality. Third, as is rigorously known for the hon- 
eycomb lattice and suggested by the LTA for the other 
lattices below, the self-dual pure Kitaev Hamiltonians lie 
in the interior of a phase rather than signifying a phase 
boundary. 

The Hamiltonians at the two Klein self dual points 
may alternatively be interpreted as possessing an en- 
larged symmetry group. The additional symmetry is gen- 
erated by the Klein duality and has Z2 characteristic. It 
thus acts in a highly nontrivial manner on spins on dif- 
ferent sites. Phases that preserve this Klein Z2 symme- 
try must contain this highly nontrivial structure; there 
is currently one known example of such a phase, the Ki- 
taev honeycomb spin liquid. If any lattice turns out to 
host a magnetically ordered phase which does not spon- 
taneously break the Klein Z2 symmetry, such a phase 
would have a complex pattern of noncoplanar spin order. 
This is unlikely, but there may also be phases which break 
the Klein Z2 symmetry but do not break too many other 
symmetries, yielding a ground state manifold that natu- 
rally splits into the two Klein Z2 broken portions. De- 
termining which or whether any of these scenarios holds 
on any particular lattice is left for future work. 

IV. EXACTLY SOLUBLE STRIPY PHASES AS 
KLEIN DUALS OF THE FERROMAGNET 

The most obvious consequence of the existence of the 
Klein duality is seen by applying the duality on the 
Heisenberg ferromagnet. At the resulting parameter 



point rj = +1, a = 1/2, the ground state manifold of 
the quantum Hamiltonian is known exactly and consists 
of simple product states, parametrized by the full SU(2) 
symmetry. The ground states may be found by taking 
a ground state of the Heisenberg ferromagnet, and ap- 
plying the rotations defined by the Klein duality on this 
magnetic order. The result is the stripy collinear mag- 
netic order. We will use the name stripy to refer to 
the FM-dual phase on lattices in any dimension, both 
to preserve the analogy to the honeycomb and also be- 
cause, as shown below, the 3D-stripy orders can have 
some "stripy" features in their own right. 

Away from the SU(2) symmetric point the symmetries 
reduce to the lattice SOC operations. The stripy ordering 
breaks the three-fold rotation symmetry, present in all 
Kitaev lattices as in Fig.jSj that simultaneously permutes 
the suitably chosen Euclidean directions a: — ^ y — >■ z — ^ x, 
the same axes on the Bloch sphere and also the Ki- 
taev bond labels a; z —>■ x. The appropri- 
ate coordinate system is set by an IrOg octahedron, in 
which the ordering is along one of the three directions 
(1, 0, 0), (0, 1, 0), (0, 0, 1) i.e. x,y,z. z-type stripy order 
has z-bonded spins aligned parallel and x or y bonded 
spins aligned antiparallel. The collinear spin axis is then 
fixed to S*^ , though the direction of the ordered moment 
will likely be determined by other effects in any material 
realization. 

The stripy orders on the various lattices share com- 
mon features but also host distinguishing characteristics. 
On the two dimensional lattices, which always appear as 
layers perpendicular to the (1, 1, 1) axis in the IrOg coor- 
dinate system, the ordering breaks the (SOC version of) 
120° lattice rotation symmetry. On the triangular lattice 
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(a) 3D-stripy order on the face centered cubic (fee) (b) 3D-stripy order on the pyrochlore lattice, 

lattice. 




(c) 3D-stripy order on the hyperkagome lattice. 



FIG. 6: 3D-Stripy orders on the fee, pyroehlore and hyperkagome lattices. Black(white) spheres represent 
up(down) spins. Bonds are colored red, green, and blue according to the Kitaev label. The ordering pattern is of 
alternating planes, here normal to z; z-type bonds (blue), i.e. those within the planes, connect spins of the same 
orientation. On the pyrochlore and hyperkagome, planes are broken up into uniformly oriented chains, with chains 
in spin-up planes oriented perpendicularly to chains on spin-down planes. On the hyperkagome, the chains are 
further broken into the linear three site clusters shown in Fig. |4] These stripy orders are exact at ri=+l, a=l/2, 

being Klein-duals of the ferromagnet. 



it is literally alternating stripes (i.e. lines of sites) of up 
spins and down spins. On the honeycomb lattice, each 
stripe is composed of the two-site clusters that lie on a 
given line; this order is also known as "IV" in the J1-J2- 
J3 literature. On the kagome lattice, stripy order gives 
the same configuration on each unit cell (is wavevector 
r) of two spins up and one spin down, meaning it is fer- 
rimagnetic with a nonzero net magnetization. At the 
exact a = 1/2 point the spins are saturated and the net 
magnetization is 1/3 that of the ferromagnet. 

In three dimensions, the 3D-stripy orders involves al- 
ternating planes of up spins and down spins. For say 



z stripy order the planes are normal to z. On the face 
centered cubic (fee) lattice, the planes are faces of the fee 
cube. On the pyrochlore lattice, the stripy order acquires 
an additional feature: spin-up planes are broken up into 
chains aligned in one particular direction, and spin-down 
planes are composed of chains aligned in the perpendic- 
ular direction. On the hyperkagome lattice this feature 
persists, and moreover the chains are broken into oriented 
linear clusters: for z-stripy order the z-type three-spin- 
chain clusters of Fig. [4] are oriented uniformly within the 
spin-up planes, and also uniformly but in a perpendicular 
orientation within the spin-down planes. The 3D-stripy 
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orders are shown in Fig. [6j 



V. LUTTINGER-TISZA APPROXIMATION 
PHASE DIAGRAMS 

Except for the Heisenberg ferromagnet and its Klein 
dual point as described above, the Kitaev-Heisenberg 
Haniiltonians are JrustrateJ^. The resulting sign prob- 
lem for quantum Monte Carlo renders their quantum 
phase diagrams, especially for the three dimensional lat- 
tices, exceedingly difficult to compute. The Hamiltoni- 
ans Eq. [l]on the various lattices are quite unique in that 
they all offer an exact solution at a nontrivial point in 
the phase diagram, the Klein dual to the ferromagnet. 
To explore the remainder of the phase diagrams we must 
use approximation methods, as we shall now describe. 

For an initial survey of the phase diagrams we employ 
the Luttinger-Tisza Approximation (LTA) , also known as 
the spherical modepSH5i]_ ^hiJe the classical version of 
a Hamiltonian has the hard constraint that the spin vec- 
tor on each site must have magnitude 5, this constraint 
may be relaxed and implemented only on average, with 
a single global Lagrange multiplier. Within this approx- 
imation the Hamiltonian Eq. [T] becomes free quadratic 
and the lowest energy configuration of the classical spins 
may simply be found by a Fourier Transform and a diag- 
onalization of the spin and sublattice indices. The LTA 
always computes a lower found to the energy of the classi- 
cal model; this inequality becomes a strict equality when 
the LTA minimum energy configuration happens to obey 
the unit length constraint. 

On a Bravais lattice and for the case with SU(2) 
spin rotation symmetry, solutions with normalized spins 
can always be constructed from the LTA minimum 
eigenvalues^. For momenta q satisfying q = —q there is a 
family of degenerate orders but even for arbitrary incom- 
mensurate momenta there are coplanar spiral solutions 
with normalized spins, with the first(second) spin compo- 
nent modulated by the real(imaginary) part of exp[i(7r]. 
However, when SU(2) spin rotation symmetry is broken 
such as by SOC, there may only be one low energy spin 
component and this approach can fail, requiring q = —q 
to construct states with unit length normalized spins. 
On lattices with multiple sites per unit cell, the LTA 
may assign different lengths to sites in the unit cell, in 
which case again we say the LTA "fails". Where the 
LTA fails, the ground state is likely to involve substan- 
tial quantum fluctuations and requires a full quantum 
analysis, though if the spins have nearly the same length 
then we expect the ordering pattern to be robust. Thus 
while the Luttinger-Tisza approximation cannot charac- 
terize non-classical phases, it is a useful first approach 
for identifying features in the phase diagram.152] 

The LTA phase diagrams are shown in Fig. [l] Here 
we discuss general features; see the appendix for details. 
Stripy phases are found surrounding the FM-dual point 
in all of the lattices; they are exact ground states at 



r]=+l, a=l/2 even within the LTA. However, the kagome 
and hyperkagome lattices exhibit an interesting frustra- 
tion: while spins are uniformly normalized at the SU(2) 
FM-dual point, away from a = 1/2 the energy is min- 
imized when spins within the unit cell are of different 
lengths. As it must by the Klein duality, this frustration 
is observed in the ferromagnet phase as well. Evidently 
for the kagome and hyperkagome, but not for the py- 
rochlore or the other lattices, even small SU(2) breaking 
within the ferromagnetic phase creates substantial frus- 
tration visible in the LTA. 

At certain points in the phase diagram, all wavevec- 
tors in the BZ offer spin configurations with the same 
minimum energy, so that the lowest band is flat. While 
subextensive degeneracies occur generically at certain 
parameter points and are expected to be completely 
lifted by boundary conditions, such extensive degenera- 
cies, marked by "Q" in Fig. [T] likely signify a new phase. 
What could the new phase be? There are only two Hamil- 
tonians hosting LTA extensive degeneracies for which the 
quantum ground state is known: the honeycomb Kitaev 
model (a = ±1) which is exactly soluble, hosting the 
Kitaev QSL with Majorana fermionic spinous; and the 
kagome Heisenberg antiferromagnet, which was recently 
found by DMRG simulations'^^ to host a QSL phase, 
consistent with a bosonic Z2 QSL^^^. The ground states 
of pyrochlore and hyperkagome Heisenberg antiferromag- 
nets, which also have LTA flat bands, are not conclusively 
known but have been proposed to be plaquette or dimer 
valence bond solids (VBS) as well as various fractional- 
ized QSLSsMZI. 

There are thus two conclusions to draw about the other 
Q points in Fig.[l] First, by the Klein duality, any lattice 
hosting a phase with no magnetic order in its AF Heisen- 
berg model also has the same type of phase surrounding 
the 77=— l,a=l/2 point, with FM Heisenberg and AF 
Kitaev exchanges. For example, the recent discovery of 
the kagome AF Heisenberg QSL then immediately yields 
the Klein dual of this QSL at the dual point; this Klein 
dual QSL will likely have distinct physical properties in 
its response to external fields. Second, by analogy with 
the known Q points mentioned above, we may guess that 
the pure Kitaev models on the hyperkagome also host a 
quantum phase with no magnetic order, either a VBS or 
a QSL. 

It is especially encouraging that the LTA flat bands 
within the honeycomb and the hyperkagome pure Ki- 
taev models arise via the same mechanism. Consider 
that LTA flat bands in the AF Heisenberg models occur 
due to the lattice specific band structure from a hopping 
model with vr flux. For the pure Kitaev models a = ±1, 
a given spin component such as hops only on z-type 
bonds. As mentioned above, for the honeycomb and hy- 
perkagome lattices, turning off y and x bonds splits the 
lattice into an extensive number of localized disconnected 
segments, as shown for the hyperkagome in Fig. m Lo- 
calization in the disjointed clusters yields the flat bands. 
Moreover, unlike for the Heisenberg case where (in the 
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relevant lattices we study) there are gapless excitations 
where the flat lowest band touches higher ones, for the 
Kitaev cases the disjointed clusters yield a band structure 
where all bands are completely flat and fully gapped, in 
the hyperkagome case also fourfold degenerate at each 
wavevector due to four clusters in the unit cell. 

Returning to the survey of the LTA phase diagrams, 
we find other regions with strong frustrations. On the 
kagome and pyrochlore lattices, over wide regions of pa- 
rameter space the LTA fails spectacularly: in the regimes 
labeled "frustrated", the unit cell in both lattices has 
two spins aligned antiparallel but with the remaining one 
(kagome) or two (pyrochlore) spins chosen to have ex- 
actly zero length by the LTA, i.e. effectively discarded. 
A related regime on the fee, a Bravais lattice, finds subex- 
tensive degeneracy involving incommensurate momenta, 
that would form spiral orders but with only one low en- 
ergy spin component cannot achieve correctly normalized 
spins across the spiral. 

Finally, on the hyperkagome the LTA finds two regimes 
with apparent magnetic order with unconventional spin 
configurations. Though in both cases the spins cross the 
unit cell are not chosen to have the same length, this is 
expected with such a large unit cell, and the LTA con- 
figurations should serve as good starting points for the 
classical and quantum Hamiltonian ground states. For 
r]—+l and a < we find that for z-type order, z-clusters 
all have the identical spin ordering "(up, down, up)", re- 
sulting in an AF state with a nonzero net magnetization 
which we thus term the "cluster-ferrimagnet" . The Klein 
dual of this order, for large a at 77=— 1, has the same "(up, 
down, up)" pattern in each 2;-cluster except clusters are 
fiipped on alternating planes, so there is zero net mag- 
netization; we term it the "cluster- AF" state. These two 
Klein dual orders are shown in Fig. [7] and Fig. [8j 



VI. SEARCHING FOR ANALOGUES OF THE 
KITAEV MAJORANA SPIN LIQUID BEYOND 
THE HONEYCOMB 

All the lattices in Figure [3] except the honeycomb have 
coordination number larger than three, spoiling the Ki- 
taev honeycomb spin liquid exact solution. However, 
similar Majorana QSL phases could still occur for the 
Kitaev Hamiltonians on the other lattices, only without 
an exact solution and with nonzero correlation length. 
Since it is generally highly difficult to determine whether 
the true ground state of a spin Hamiltonian forms a QSL, 
we will not attempt to answer this question. Instead, we 
will study possibilities for similar Majorana QSLs on the 
other lattices using an appropriate choice of mean field. 

The exact solution of the Kitaev honeycomb model in 
terms of Majorana fermion operators is a specific case 
of a Schwinger fermion decomposition mean field, which 
becomes exact for this modeP^"^. To search for similar 
Majorana QSLs on the other lattices, we thus employ 
this mean field. Spins are decomposed into bilinears in 




FIG. 7: Magentic order of the cluster-ferrimagnet 
state on the hyperkagome. Black(white) spheres are 
up(down) spins; bonds are colored according to Kitaev 
label. Here shown for z ordering, notice that the z-type 
clusters (blue), lying in planes normal to 2, all have 
"(up, down, up)" spin configurations. This 
configuration has nonzero net magnetization. Note that 
the cluster-ferrimagnet order is Klein-dual to the 
cluster-AF order. 




FIG. 8: Magentic order of the cluster-AF state on 
the hyperkagome. Black(white) spheres are up(down) 
spins; bonds are colored according to Kitaev label. Here 
shown for z ordering, notice that the z-type clusters 
(blue), lying in planes normal to z, have "(up, down, 
up)" spin configurations on even-numbered planes, and 
the opposite "(down, up, down)" spin configurations on 
odd-numbered planes. This configuration has zero net 
magnetization. Note that the cluster-AF order is 
Klein-dual to the cluster-ferrimagnet order. 

four Majorana species ;^o^i'2,3 

^''^zxV with {x^/l^-^"''- (14) 

This mapping is exact under the single fermion occu- 
pancy constraint XiXlXiXi = 1/4- On the honeycomb 
lattice this constraint commutes with the pure Kitaev 
Hamiltonian, but that does not occur on the other lat- 
tices. The Z2 gauge freedom in defining the Majorana 
operators enables a choice of attaching gauge transfor- 
mations to the physical symmetry operations, called a 
projective symmetry group (PSG)|Sil see the Appendix 
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for details. The PSG of the Kitaev honeycomb model 
was previously studiecP^ and determined to be flux-free, 
with (xi, X2j Xa) transforming as a pseudovector and each 
bond permitting Majorana bilinear expectation values 
only for two Majoranas of the same species a, yielding 
a total of three mean field parameters: 

U° = Wo, K = Ua, u''/" = Ub (15) 

with 

^ (16) 

where 6 is a bond. Uf, is set to zero for the pure Kitaev 
model. The resulting mean field Hamiltonian is 

Hmf = -hT,i,v[^.,aSign[v]v^^^^ixfxf 
^ 7] {{1 - \a\) - 2aS^) uo 
^" = Ea ■^7^7 = ?? ((2 - 2a - 2\a\)ub ~ 2aUa) (17) 

where J° is the coupling of spin component a on a 7- 
bond, i is a site, v[i] are the bonds of site i and sign[u] is 
the orientation of the bond v within the PSG. This orien- 
tation determines how operators on the bond transform 
under symmetries. On the honeycomb, bonds are ori- 
ented from sublattice A to B. The bond orientations used 
in the PSGs for the triangular and the kagome lattice are 
depicted in Ref . [531 (though that work dealt with bosonic 
QSLs, the bond orientation diagrams are the same). For 
the triangle, it is known as the zero flux PSG. For the 
kagome, this zero flux PSG is known as y/3 x \/3 or 
Qi — —Q2- The PSG analysis for this type of mean 
field has not been successfully carried out on the 3D lat- 
tices; the pyrochlore does not appear to give a unique 
decompositiorP2l. On the hyperkagome however one of 
the four spins in each tetrahedra is removed, so we can 
consistently choose the orientation A B C —?' A 
within a triangular face in Fig. |6c] giving a unique PSG. 

The mean field Hamiltonian Hmf is a free Majorana 
bilinear Hamiltonian, so its ground state is immediately 
known by computing its band structure. The qualitative 
properties of this band structure carry the primary infor- 
mation, though the band structure energy scales contain 
the unknown mean field parameters u. The parameters u 
can be determined self-consistently from the band struc- 
ture by computing the Majorana propagator, as a Mat- 
subara frequency integral of the inverse of the frequency- 
dependent Hamiltonian kernel; computing it analytically 
is nontrivial on a lattice with more than two sites per 
unit cell. We have carried out the self consistency com- 
putation on the triangular and honeycomb lattice, using 
the Kitaev-type majorana flux-free PSG which is defined 
for both, and find that the mean fields evolve with a 
smoothly away from the Kitaev limit, with no first order 
transitions. 

Regardless of the exact values of the mean field pa- 
rameters, choosing the mean field to be analogous to the 
Kitaev honeycomb QSL already determines key proper- 
ties of the resulting states on the various lattices. First, 



all the lattices except for the honeycomb possess cycles 
with an odd number of bonds, such as triangles; this im- 
mediately requires the Kitaev Majorana mean field to 
spontaneously break time reversal symmetrji^. These 
time reversal broken spin liquids might not display typ- 
ical characteristics of time-reversal broken states. For 
example, on the triangular lattice, even though time re- 
versal as well as 2n/6 rotation each independently flip 
the flux pattern in triangular faces, the combined opera- 
tion of time reversal with 27r/6 rotation is still preserved 
as a single symmetry operation, so the Hall conductance 
vanishes. Second, lattices with an odd number of sites 
per unit cell necessarily have a spinon Fermi surface; 
the even-unit-cell lattices of pyrochlore and hyperkagome 
may or may not host gapped spinous. 

Third, certain qualitative features of the band struc- 
ture are determined by the choice of mean fields, such as 
the consideration of only nearest neighbor bonds and the 
PSG. There are four Majorana fermion species per site; 
for a pure Kitaev Hamiltonian, x^''^'^ have bands related 
to each other by the 120° SOC combined spin-spatial ro- 
tation, while x° has a generally different dispersion. For 
the honeycomb model, x° has a Majorana analogue of the 
Dirac cone, i.e. relativistic with zero mass, while x^'^'^ all 
have completely flat bands separated from zero energy by 
a complete gap. The kagome lattice x^'"^'^ also has a flat 
band but it lies at zero energy i.e. at the Fermi energy, 
yielding the Fermi surface which necessarily arises here. 
The flatness results from a localized unpaired Majorana 
mode on one of the three sites in each unit cell; but since 
the remaining two sites form a line spanning the lattice, 
they disperse and the other bands are not flat, touch- 
ing zero energy along lines in a quasi-lD spectrum. For 
the pyrochlore even qualitative statements cannot be cur- 
rently made, since as mentioned above, there is no special 
choice of minimal flux PSG. On the hyperkagome with 
bond orientations as described above, x^ has some gap- 
less subextensively degenerate modes (such as from F to 
M); but x^'^'"^, hke for the honeycomb, have completely 
flat bands. These arise, as previously mentioned, because 
both the honeycomb and the hyperkagome fragment into 
extensively many disconnected clusters when only bonds 
of a single Kitaev label are kept. However, while the hon- 
eycomb clusters have an even number of sites and hence 
can form two fully gapped bands, separated from zero 
energy, the hyperkagome clusters have an odd (three) 
number of sites; each cluster always has one energy band 
at zero energy and hence x^'^'^ are gapless. 



VII. OUTLOOK 

On the honeycomb lattice, the roles of the SU(2)- 
symmetric Heisenberg coupling and the SOC Kitaev cou- 
pling are distinct and clear: Heisenberg exchange yields 
magnetic order, Kitaev exchange yields the exactly solv- 
able QSL phase. The natural interpolation between the 
two limits, that would occur if the couplings arise in irid- 
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ium oxide compounds, is consistent with this framework: 
the intermediate region simply holds more magnetic or- 
der. However, as we have discussed above, generaliza- 
tions of the Kitaev coupling naturally arise in iridium 
structures and other geometries of edge-sharing octahe- 
dra on many other lattices, motivating the study of the 
phase diagrams of Eq. [T] on these various 2D and 3D 
lattices. Beyond the honeycomb, the roles of the two 
exchanges begin to break down. 

The effect of lattice geometry on the "frustration" of a 
lattice is quite different for the two terms; the Hamilto- 
nian and the lattice determine the frustration together, 
not independently. More surprisingly, even in cases when 
the Heisenberg Hamiltonian appears highly frustrating, 
interpolating between the AF Heisenberg and the Kitaev 
limits, we find a phase which occurs on all the lattices 
and which is exact and fluctuation free at a certain pa- 
rameter point. Subtle interplays of different magnetic 
couplings, rather than a monotonic "frustration" mea- 
sure, seem to be at play. The intermediate stripy phase 
is exact by virtue of being related to the ferromagnet, 
through a duality that emerges through the SOC on the 
t2g orbitals microscopically generating the Hamiltonian. 

The Klein group structure of the mapping between 
Hamiltonians (a duality) is in some sense highly spe- 
cific to these quantum chemistry considerations but in 
another sense, as a mathematical object Z2 x Z2, quite 
universal. The duality transformation it generates is in- 
teresting here for another reason: while most dualities fix 
a single self-dual Hamiltonian and map the two regimes 
on either side of that point, with qualitatively different 
features, to each other, the Klein duality is different. It 
admits two self-dual Hamiltonians, which seem to gen- 
erally lie in the interior of a phase. And it acts in a 
complicated way on spin and spatial indices, making its 
action as a Z2 symmetry operation highly nontrivial. 

Regarding possible experimental significance of Hamil- 
tonians arising from strong SOC, it is important to ob- 
serve that the Kitaev couplings naturally occur in a man- 
ner more subtle and constraining than naive symmetry 
considerations would suggest: for example, Kitaev inter- 
actions can arise for iridium ions on the fee but not on the 
simple cubic. Computations of the quantum phase dia- 
grams on the various lattices, especially the pyrochlore 
and hyperkagome, will pave the way towards predictions 
and comparisons with experimental results. 
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Appendix A: Details on the Luttinger-Tisza 
Approximation phase diagrams for the various 
lattices 

Here we give the results of the Luttinger-Tisza Ap- 
proximation (LTA) analyses in more detail. The classi- 
cal Kitaev-Heisenberg model on the triangular lattice has 
been recently studied"**^, and the phase diagram for the 
honeycomb is already knowrPl we begin by reviewing the 
results for the triangular and honeycomb lattices for the 
sake of direct comparison with the other lattices. Note 
that on all the lattices, the pure Kitaev points a = ±1 
always host degeneracies in the BZ; except where noted 
below, these are just subextensive degeneracies, which we 
expect will be lifted by small perturbations or by bound- 
ary conditions for a finite system. 

Triangular lattice: AFM Heisenberg, rj = +1. The 
pure Heisenberg antiferromagnet a = hosts ordering at 
wavevector K (BZ corner), corresponding to 120° order 
with a tripled unit cell. This ordering can also be seen 
by explicitly working with the enlarged unit cell. For 
< a < 1/5 the ordering wavevector is incommensurate 
as it migrates from K to M (BZ edge midpoint). For 
— 1 < a < the ordering wavevector is also incommen- 
surate, moving from K in the direction of F. Only one 
spin component has minimum energy for a given momen- 
tum here, preventing the construction of any order with 
unit length spins such as the usual coplanar spiral. The 
antiferromagnetic Heisenberg limit 120° order actually 
has point topological defects (vortices), with Z2 charac- 
ter; the recent triangular lattice classical model study in- 
cluded classical Monte Carlo computation which found 
they may be stabilized into a vortex lattice in this in- 
commensurate ordering phase. Finally, for 1/5 < a < 1, 
we find stripy order at wavevector M. FM Heisenberg, 
77 = —1. For — 1 < a < 3/7 there is ferromagnetic order. 
For 3/7 < a < 1 continues the phase of incommensurate 
ordering wavevector between F and K from 77 = -1-1. 

Honeycomb lattice: For pure Kitaev exchange a — 
±1 there are extensive degeneracies. AFM Heisenberg, 
rj — +1. For —1 < a < 1/3, simple Neel order (wavevec- 
tor F) is stabilized, as the honeycomb is bipartite. For 
1/3 < a < 1 there is stripy order at M. FM Heisenberg, 
rj — ~1. Because we have a quadratic spin Hamiltonian 
and a bipartite lattice, this regime can be mapped to 
?7 = -|-1 by flipping all spin components on one of the two 
sublattices. For — 1 < a < 1/3 there is a ferromagnet, 
for 1/3 < a < 1 there is zigzag order at M. 

Kagome lattice. AFM Heisenberg, rj = +1. The 
a — Heisenberg antiferromagnet hosts extensive degen- 
eracy. For < a < 1 there is stripy order (wavevector 
F), here ferrimagnetic. But while spins are uniformly 
normalized at the SU(2) point a = 1/2, away from 
a = 1/2 the energy is minimized when spins within the 
unit cell are of different lengths, suggesting the order- 
ing becomes increasingly frustrated. For example near 
a — 1/2, if two of the spins within the unit cell are cho- 
sen to be normalized with magnitude -1-1, the third has 
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magnitude 2(-i+a) > or roughly — 2+2a. l<or 

— 1 < a < the LT method fails: two spins within a tri- 
angle plaquette align antiparallel along the axis of their 
bond label 7, but the third is so frustrated that its mag- 
nitude is set to exactly zero, also giving spurious subex- 
tensive line degeneracies from T to M. FM Heisenberg, 
77 = —1. a = 0.5, the Klein rotation of the Heisenberg 
antiferromagnet, again hosts extensive degeneracy. The 
ferromagnet exists over a wide range — l<a<l/2 but 
with the same strong frustration as for the stripy phase: 
away from the Heisenberg limit a = 0, the energy is 
minimized by having spins of different lengths within the 
unit cell. For example, at a = —2/5, the unit cell has 
all spins aligned but with one spin at 2/3 the magni- 
tude of the others. The regime 0.5 < a < 1 continues 
the antiferromagnctic-Kitaev regime of 77 = +1, with two 
oppositely aligned spins and the third spuriously set to 
zero. 

Face centered cubic (fee) lattice. AFM Heisen- 
berg, T] = +1. 0<a<l hosts the 3D-stripy order, 
antiferromagnetic along a cubic axis, which for the fee 
BZ is at wavevector X (center of one of three square 
faces of the BZ). For order, stripy-ordered spins align 
along S'^.s For — 1 < a < 0, there is a line degeneracy 
from XtoW (i.e. along diagonals of a square face of the 
BZ). The square face lies at some Euclidean direction, 
say z. The two perpendicular spin components (here S'^ 
and S^) participate at X, with one of the two stabilized 
as the wavevector moves towards a corner W. The set 
of minimal energy spin configurations, with normalized 
spins, therefore contains the wavevector X stripy order 
though with spins imiformly rotated to the , plane. 
The orders at generic degenerate g = (1 — r)X + rW 
would usually form a spiral but here only have mini- 
mized energy for one of the three spin components, so 
they cannot stabilize any order with unit length spins. 
For the constrained classical Hamiltonian and the quan- 
tum Hamiltonian, the higher energy spin component (s) 
are likely to be mixed into the ground state, so it ap- 
pears the LT method hints at spiral order in this regime. 
FM Heisenberg, rj = —1. The ferromagnet is stable for 
-1 < a < 1/2. The region 1/2 < a < 1 is part of the 
spiral-like phase of 77 = +1, forming a single phase stabi- 
lized by antiferromagnetic Kitaev interactions regardless 
of the sign or presence of Heisenberg exchange. 

Pyrochlore lattice. AFM Heisenberg, 77 = -1-1. 
The a = Heisenberg antiferromagnet hosts extensive 
degeneracy. Stripy order (wavevector F) is stable for 
< a < 1. For -1 < a < the LT method fails: 
the preferred states have two of the four spins in the unit 
cell pointing antiparallel to each other, with the other two 
spuriously set to zero. The (subextensive) degeneracy ex- 
tends across planes in the BZ, specifically the high sym- 
metry planes containing the gamma point and the centers 
of (both hexagonal and square) faces. FM Heisenberg, 
T] = —1. a = 1/2, the Klein rotation of the Heisen- 
berg antiferromagnet, again hosts extensive degeneracy. 
The ferromagnet is stable for —1 < a < 1/2. Strong 



antiferromagnetic Kitaev exchange atl/2<Q;<l con- 
tinues the 77 = +1,Q! < regime of plane degeneracies 
and two missing spins. There are subextensive degenera- 
cies at a = ±1. Note that unlike the kagome case, here 
even away from SU(2) points, the ferromagnet and stripy 
regimes have all spins of unit length. 

Hyperkagome lattice: AFM Heisenberg, r] = +1. 
The a = Heisenberg antiferromagnet hosts extensive 
degeneracy, as do the two pure Kitaev limits a = ±1. 
Stripy order (wavevector F) nominally exists for < 
a < 1, but with the same strong frustrations (except 
for the SU(2) point a = 1/2) as for the kagome. For 
say z ordering, the four spins in the middle of their re- 
spective 2;-type clusters are chosen smaller(larger) than 
unity for a < 1/2(q! > 1/2). The region with anti- 
ferromagnetic Kitaev exchange — 1 < a < hosts a 
magnetically ordered phase which we term the cluster- 
ferrimagnet. In the 2- type cluster-ferrimagnet, ^-clusters 
all have the identical spin ordering "(up, down, up)", 
resulting in a state with a nonzero net magnetization. 
The LTA gives the spins in the middle of each cluster 
a somewhat larger, unphysical, normalization, but this 
magnetic pattern is likely to be robust. FM Heisenberg, 
77 = — 1. a = l/2, the Klein rotation of the Heisenberg 
antiferromagnet, hosts extensive degeneracy like in the 
pyrochlore. (As already mentioned, the a = ±1 points 
do as well.) The ferromagnet over — 1<q;<1/2 again 
has different length spins within the unit cell, except for 
the SU(2) Heisenberg point a = 0. For 0.5 < a < 1 
there is the cluster- antiferromagnet, which is the Klein 
dual of the cluster-ferrimagnet. The two orders appear 
to be separated by the likely quantum phase, associated 
with extensive degeneracy, at 77a = -|-1. In this order, like 
for the cluster-ferrimagnet, the middle spin in each clus- 
ter points opposite to the other two and prefers to have 
a larger magnitude. Unlike the cluster-ferrimagnet, the 
cluster-antiferromagnet has no net magnetization since a 
pair of clusters is flipped relative to the other pair, as can 
be seen by the Klein duality. 



Appendix B: Mean fields and the projective 
symmetry group (PSG) 

The choice of partons (here Majorana fermions) with 

which to decompose the spin is not unique; there is ac- 
tually full gauge SU(2) freedom in defining the parton 
operators. When parton mean field hopping and pair- 
ing terms gain nonzero expectation values, the freedom 
is partially lost so that the mean field Hamiltonian is 
only invariant under a global Z2 transformation. This 
is known as the invariant gauge group (IGG). Beyond 
the mean field level, the partons couple to a dynamical 
gauge field, with the gauge group of the IGG. Hence when 
IGG= Z2, the spin liquid is known as a Z2 spin liquid. 

The mean field ansatz may be chosen to preserve as 
many of the symmetries of the original Hamiltonian as 
possible, say a set of symmetries S. While the mean field 
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ansatz (the collection of mean field amplitudes on the 
lattice) may not be immediately invariant under a sym- 
metry in S, we may choose to always attach an SU(2) 
gauge transformation to the symmetry, such that the 
mean field ansatz is invariant under their combination. 
Relations of S (operators that multiply to the identity) 
must have associated gauge transformations that multi- 



ply to an element of the IGG, ie a global 1 or —1 sign in 
the case where IGG= Z2. This collection of attached real 
and gauge symmetry operations that leave the mean field 
Hamiltonian invariant is known as the projective symme- 
try group (PSG). It is a necessary part of the description 
of any mean field ansatz, and distinguishes different spin 
liquid phases even when they have the same gauge group. 
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